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The purpose of this article is to address the dynamics of an interacting Bose-Einstein conden- 
sate confined in coupled one-dimensional Landau-Zener arrays under the influence of disorder and 
harmonic confinement. In particular, we concentrate in studying the interplay of disorder and in- 
terparticle interaction on the transfer of atoms depending on the speed of Landau-Zener sweeps. A 
dynamical phase diagram summarizing the final situation across ground state and inverse sweeps is 
given in terms of the effect of disorder, interaction and the speed of the sweeps. 



The most important agents that determine the dynam- 
ical behavior, and thus the properties, of a quantum sys- 
tem are their structure and the interactions among their 
constituents. Regarding the influence of the confinement 
potential, the seminal work by Anderson [1 investigating 
the conductance of an electronic system at T = showed 
the striking effect of a disordered structure on the change 
in the conduction in a lattice. On the other hand, the 
understanding of bosons in disordered media provided by 
the work of Fisher [2] through the Bose Hubbard model 
demonstrated the existence of superfluid and Mott insu- 
lating phases as a function of the ratio between interac- 
tions and hopping strengths, restricting in the Mott insu- 
lating case the transport across the media. The structure 
itself works as an energetic landscape where the atoms 
or particles in general, are allowed to move under the in- 
fluence of the interactions. When particles are trapped 
during their evolution into long-lived inhomogeneous sta- 
ble or met ast able states, a Mott insulating or glass phase 
can emerge respectively [3]. Apart from the conditions 
that define the accessible quantum states consistent with 
the closed system, there are additional factors that may 
prevent the particles to explore the available quantum 
Hilbert space in their evolution towards the stationary 
state. We refer to dynamical restrictions imposed exter- 
nally during their time evolution. Such protocols are the 
so called quantum quenches, and can effectively change 
either the interactions among the system constituents 
[4] or the confinement potential. Landau-Zener (LZ) 
phenomena [5] are ideal scenarios where dynamical con- 
strains can be recreated since adiabatically or at a given 
rate, the energy landscape or potential energy surface is 
continuously changed, enhancing the particles to modify 
its natural behavior. 

Even though the intrinsic inhomogeneity created by 
the harmonic confinement in the experimental realization 
of ultracold atomic gases, most of the condensed matter 
and solid state phenomena can be reproduced in the lab- 
oratory by adapting the external fields where the atoms 
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FIG. 1: (Color online) Schematic representation of a quan- 
tum Bose liquid in chain of double-well traps along a sweep 
from to A (figures (a) and (b)). Disorder is introduced by 
changing randomly the potential depth at each site a quantity 
that varies between what we called strong disorder and weak 
disorder (dotted and solid lines in figure (c)). 



move and by adjusting the interactions among the parti- 
cles in diluted conditions. Particularly, when the atoms 
are confined in the so called optical lattices, the insulating 
and glass phases mentioned above can be materialized by 
modifying the local interactions and setting up a spacial 
disorder respectively. Bosons [6HTU]. fermions [11 and 
mixtures of them p~2j [13] confined in random potentials 
have demonstrated to exhibit such an insulating and glass 
phases, being dimensionality another factor participating 
in the appearance of localization phenomena [TJJ [15] . 

The distinctive characteristics that establish the na- 
ture of which a quantum phase is present in a Bose fluid 
at zero temperature are measured in terms of their su- 
perfluid fraction, the existence of a gap in the excitation 
spectra, the behavior of the compressibility, the density 
and the fixed (integer or semi-integer) number of bosons 
at each lattice site. In this work we study the many 
body LZ generalization proposed by Y.-A. Chen [16] et 



at. consisting of a Bose Einstein condensate confined 
in two-coupled chains in ID where the potential depths 
defining the lattice sites are linearly modified in time to 
sweep from an initial energy difference —A to a final dif- 
ference A. In addition to such a time dependent poten- 
tial, we incorporate the presence of disorder by changing 
the potential depth at each lattice site. 

First we establish the model that describes the system 
under consideration. The many body LZ generalization 
is represented through the Hamiltonian [T6HT8] : 

i i 
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Hi = -J 1 - (b\ R b i L + h.c.) + A (hi, R - hi, L ) 

H v i= \ E KMXA»> (!) 

v=R,L 

where b\ v and b i v are the usual creation and destruc- 
tion Bose operators with its associated number operator 
hi )V . J 'I and J x are the intra- and inter-chain coupling 
energies, and U is the onsite inter-particle interaction 
strength. While labels i and j designate sites, for v\ R 
and L are chosen to denote right and left chains respec- 
tively. The presence of disorder is introduced by changing 
in a random manner the potential depth at each lattice 
site. Such a variation is effectively incorporated in the 
model by scaling the number of particles with a random 
number (white noise) 5i iV in an interval [—5, 5] from with 
mean 0. Therefore, the Hamiltonian describing disorder 
at each site has the form, 

H?= J2 5 ^<- ( 2 ) 

v=R,L 

Similar disorder can be realized experimentally by an 
speckle potential as in p~9j [20] . In our analysis we 
have considered that disorder is the same in both tubes 
Si = 5i,R = Si,L- The effect of the harmonic confinement 
associated to the attainment of the condensate is outlined 
by means of a Hamiltonian term: 

Uf = ^2 V ( x i ~ x o) 2r hi,v, (3) 

u=R,L 

where V modulates the curvature of the harmonic con- 
finement and xo denotes the center of the trap, xi = 
(2xq /N x )i. Thus, the full quantum model considering 
the effect of the trap and disorder in the LZ scenario is: 

h = n LZ + ( n ? + n i) - ( 4 ) 

i 

Numerical calculations for realistic number of particles 
become impracticable since Hilbert space scales as g 2Nx 
and thus dynamical evolution for Hamiltonian Q re- 
quires to diagonalize a matrix of size g ANx , where N x is 
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FIG. 2: (Color online) The density of transferred atoms 
to the right side of the chain pR after a ground state 
sweep as a function of the speed of the sweep a and 
the position along the system x, in dimensionless units. 
Parameters are: U — 10.0, V — 15.0/xq, xq — 50; 5 = 
0(a, no disorder), 4. 0(b, weak disorder), 20. 0(b, strong disorder), 
see main text for details. 

the number of sites and g is the filling factor. As shown 
in [16] . the number of sites is ~100. Then, to simulate 
the system with both disorder and the effect of the har- 
monic confinement, we use the decoupling approximation 
per site [18], such that the Hamiltonian becomes a sum 
of local contributions per site Hf z . 

i 

where the LZ Hamiltonian is given by, 

u=R,L 

(6) 

being ^r/l order parameters determined self- 
consistently, such that they minimize the ground 
state energy of the system for a given set of all the other 
parameters. The macroscopic state of the system can be 
written as: 

|*>«|&)®---®|^.> = ®&|&), (7) 

where \<pi) is a local many-particle state. The time evo- 
lution of the many-particle macroscopic state is given 
approximately by: 

|*(t„+i)> « e- in ^ st \^(t n )} = ^e-^^lMtn)) 

(8) 

The time evolution is done by first finding self 
consistently the order parameters that minimize the 
ground state for all other parameters fixed, where self- 
consistently the Hamiltonian depends on time linearly 
via A(t) = at. We do this for each value of a, and we 
proceed with the same calculation as in [18 , for the 
time evolution at each site i. We should mention that 
in our calculations all the parameters are referred to the 
tunnelling coupling parameter J = jH = J x . Here and 
henceforth we use the notation U = U/J, V = V/J, 
A = A/J and 5 = 5/ J. One of the observables that we 
analyze is the local density of particles in the right tube 
after a LZ sweep, p R (i,rj,6) = (^i(rjJ)\n iiR \4>i(rj,6))/g, 
for a given set r] = {a, £/, V} and disorder amplitude S 
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FIG. 3: (Color online) The density of transferred atoms to 
the right side of the chain pR (a-c) and the compressibility 
kr after a ground state sweep as a function of the position 
along the system x, for different disorder strengths shown in 
the bar on the right, in dimensionless units. Parameters are: 
27r/a = 10.0 (adiabatic transfer), V — 15.0/xq, x o = 50, 
U = 0.1(a,d), 10.0(b,e), 20.0(c,f). 

where \4>i(r), S)) is the final state after the LZ sweep. We 
have considered a filling factor of g = 2, similar to the 
experimental conditions of [16], which means up to two 
particles per link in the tubes at a given position Xi. In 
order to have meaningful quantities in terms of the dis- 
order amplitude 5 : we consider sets of 40 realizations of 
random numbers for a given disorder amplitude, generat- 
ing the disorder contribution to the Hamiltonian ([2|. Af- 
ter the time evolution is done, we average the realizations 
for each disorder amplitude generated of the observables 
p and ft for each site i. A collection of some of these re- 
sults are shown in Fig. |2][4j To support our findings, we 
go beyond, performing calculations with sets of random 
numbers with an order of magnitude larger of realizations 
(~ 200) and we found qualitative agreement with smaller 
sample sizes. We consider disorder amplitudes in the 
interval 2 (weak disorder) < 5 < 20 (strong disorder). 
The LZ sweeps have been made from ti = —Ao/a to 
tf = Ao/a for the ground state sweeps and vice- versa 
for the inverse sweeps with Ao = 20. In our simula- 
tions we considered N ~ 100 atoms in the system, with 
N x = 64. 

As shown in Fig [2] the final density after the sweep in- 
herits the effect of the trapping potential, similar to an 
inverted parabola depending on the strength V. As one 
increases V, less atoms can be contained in the trap. We 
also consider the final state to determine the quantum 
fluctuations in the particle number per site and the local 
compressibility, K R (i,r),5) = f3($i(ijJ)\hl R \0i(ri,5)) - 

(4>i(r],5)\hi iR \4)i(r],S)) 2 ), with f3 = k B T. In overall, we 
found that as we increase the disorder amplitude 5 along 
the trap, the system becomes more homogenous. This 
comes as no surprise since disorder acts in the diagonal 
terms in the Hamiltonian, therefore as we increase the 
disorder amplitude, the atoms see on average less curva- 
ture of the harmonic confinement. Thus the system size 



increases, and the density flattens, see Figj2] and Fig(3] 
panels (a-c). In Figj3]we show the effect of the trapping 
potential for the state which is approaching the limit 
of adiabatic transfer (a — >• 0) for ground state sweeps 
(2n/a = 10) as a function of the disorder amplitude for 
(a,d) weak interaction [U = 0.1], (b,e) medium inter- 
action [U = 10.0], (c,f) strong interaction [U = 20.0]. 
Consistent with has been analyzed in [T6UT8] we found 
the breakdown of the regular LZ scenario [5] where the 
maximum transfer occurs in the fast sweep limit. We 
find consistent behavior with the LZ scenario for inverse 
sweeps. For weak interaction, disorder destroys the Mott- 
Insulating (MI) phase occurring in the center of the trap 
for 5 = 0, see panels (a) for the density [21] and (d) for 
the compressibility [22], here quantum fluctuations start 
to increase with disorder moving the region of the cen- 
ter of the system away from approximate integer filling. 
However, as we increase disorder the maximum transfer 
decreases consistent with localization driven by the An- 
derson mechanism [1, 19 , 20 . In panel (d) we can see how 
the compressibility compensates the effect in the density, 
due to the curvature away from the center of the trap for 
weak disorder, this compensation mechanism vanishes as 
the MI region disappears. As we increase interaction in 
the adiabatic limit from its medium value (b,e) to the 
strongly interacting limit (c,f) we find that the effect of 
disorder starts to saturate and the effect of the interac- 
tion prevents the disorder from affecting the system and 
as one would expect, the compressibility over all the sys- 
tem decreases. Provided that the interaction is below 
some threshold the system is affected by disorder and 
the maximum transfer of particles is controlled by the 
Anderson mechanism in the adiabatic limit, suppressing 
the transfer locally. However, as we move away from the 
adiabatic limit, the behavior of the system becomes more 
interesting. 

Another relevant set of global measures are the spa- 
tial averages, ^(rj,S) — Y2i€(hVifi)/N x , where £ can be 
either the compressibility n or the density p. Using the 
above definitions, one can construct the following quan- 
tity which measures the net effect of disorder on average 
along the system, 7^ = £(77, 5 7^ 0)/£(r/, 5 = 0) — 1. 

Due to the global effect of disorder along the trap, a 
good way of accounting for the effect of disorder is via 
r ) PR and j KR , which depending on the speed of the sweep 
a, present non-monotonic character and are strongly in- 
fluenced by the many-body interaction, see Fig. [4] In the 
weakly interacting regime, Figj4](a,c), we find that in the 
diabatic limit (27r/a — » 0), fast sweeps, the disorder ini- 
tially suppress both the population transfer and the fluc- 
tuations up to 30%. As we slow down the sweeps a max- 
imum occurs (an enhancement of 20%) in the population 
of right tube, while the effect in the compressibility van- 
ishes. Slowing down even further returns the system to 
have on average a small suppression effect on the density 
while the mean compressibility is strongly suppressed. 
Increasing the interaction, Fig. |4](b,e), leads to a point 
where the peak in the effect of the density sharpens and 
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FIG. 4: (Color online) The average effect along the trap of 
disorder in the density of transferred atoms to the right side 
of the chain ^ PR (a-c) and the effect on the compressibility 
7^ r after a ground state sweep as a function of the speed 
of the sweep a, for different disorder strengths shown in the 
bar on the right, in dimensionless units. Parameters are:y = 
15.0/xg, x = 50, U = 0.1(a,d), 10.0(b,e), 20.0(c,f). 



enhances the average density pR up to 40% for 27r/a ~ 1 
[here the maximum in pR ~ 0.5, see Fig|3](b) ], while the 
effect in the fluctuations with respect to the absence of 
disorder is marginal. Afterwards, as the effect in the den- 
sity becomes marginal the fluctuations are suppressed, as 
in the case of weak interaction. On the other hand fluctu- 
ations are strongly suppressed in the adiabatic limit and 
strongly enhanced for fast sweeps Fig. [4](e). However for 
strong interaction, Fig. [4](c,f), the scenario completely 
changes and the effect of disorder for fast sweeps is to ini- 
tially enhance the transfer (One must note that although 
the enhancement is large the population transfer is ex- 
ponentially small). Decreasing the speed suppress the 
density and then as we approach the adiabatic limit the 
peak in the effect on the density moves to slower speeds 
until it disappears. The fluctuations, Figj4](f) in general 
are being suppressed, in agreement with Fig. [3] The 
enhancement or suppression effect in the density with re- 
spect to the speed a has to do with the fact that globally 
disorder makes more homogenous the system, see Fig. [2] 
which corresponds to 5 = 0.0, 4.0, 20.0 and U = 10.0 and 
(b,c) Fig. [4j 

In the case of inverse sweeps the behavior consider- 
ing disorder is the opposite with respect to 2ir/a in the 
ground state sweeps, meaning that for fast sweeps the 
suppression effect of disorder in the density is stronger 
for weak interaction, and there exists also a saturation 
behavior for strong interaction in the maximum transfer 
(One must note that for inverse sweeps the population 
transfer in the adiabatic limit is exponentially small). 
In essence, the adiabatic limit behavior for ground state 
sweeps and the fast sweeps limit behavior for inverse 
sweeps are symmetric with respect to the effect of dis- 



order. We summarize all of our findings in a qualitative 
dynamical phase diagram in terms of the effect on the 
mean transfer density, see Fig. [5] The net effect of disor- 
der and interaction in the mean density, can be explained 




FIG. 5: (Color online) Qualitative dynamical phase diagram 
of the system for either ground state or inverse sweeps depend- 
ing on the interaction strength U and the disorder amplitude 
5, units are dimensionless. The boundary between disorder 
dominated transfer (DDT) and interaction dominated trans- 
fer (IDT) is given by the black lines. The behavior concern- 
ing the speed of the sweeps is given by the black dotted line, 
the boundary changes from (a) 2n/a ^> 1 [adiabatic limit] 
to (b) 27r/a ~ 1 to (c) 27r/a <C 1 [fast sweeps]. The gray 
(pink) areas denote the region where the effect on the den- 
sity is marginal, the brighter area (yellow) denotes the region 
where enhancement occurs via disorder and the darker areas 
denote suppression either via disorder (lower right) or via in- 
teraction (upper left). The white dots denote the approximate 
boundary between weakly interacting ( WI) , strongly interact- 
ing (SI), weak disorder (WD) and strong disorder (SD). The 
white dashed line denotes the critical value of U c for which 
interaction dominates and saturates the maximum possible 
transfer independent of disorder, see Fig. [3](c). 



in terms of suppression via interaction effects [upper dark 
region (blue) in Fig. [5] where a saturation effect with re- 
spect to disorder occurs. While disorder depending on its 
strength 5 and the speed of the sweep can have an aver- 
age suppression effect [lower dark (blue) region in Fig. [5] 
or average enhancement effect [middle bright (yellow) re- 
gion, in Fig. [5] . The local suppression effect [see Fig. [3] 
via disorder is consistent with the Anderson localization 
mechanism and strongly affects the dynamics, as well as, 
the effective confinement seen by the condensates of the 
tubes. Interference effects homogenize the system in the 
dynamical evolution leading to the picture of "glassy dy- 
namics" . 
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